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Abstract 

The resonant and non-resonant amplification of electromagnetic inhomogeneities is investi- 
gated within a phase of coherent dilaton oscillations. We show, in some explicit toy model, 
how the critical energy density bound applied to the amplified fluctuations of the massless 
gauge fields can constrain the dilaton mass. 
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1 Introduction 



The coherent oscillations of a massive scalar field minimally coupled to gravity are a possible 
source of many cosmological problems which can result in a quite stringent set of bounds 
on the value of the scalar mass. The implications of the so-called "Polonyi potential" were 
analysed in different contexts |l|, 0, 0, ^. 

The energy density stored in the coherent oscillations red-shifts faster than the surround- 
ing relativistic plasma and its contribution might become soon quite significant depending 
upon the initial amplitude of the scalar field at the beginning of the oscillating phase. If 
the initial amplitude is of the order of the Planck mass (as it seems natural to require for 
interactions of gravitational strength) the red-shifted energy density of the oscillating field 
might have catastrophic consequences for the big-bang nucleosynthesis. 

At the nucleosynthesis curvature scale ifiv(tiv) the background temperature should be 
T/v ~ 0{lMev), implying a bound for the reheating temperature associated with the scalar 
field decay and a lower bound for the scalar field mass which should be heavier than lO^Gef . 
The decay of the coherent oscillations might have dangerous consequences even before the 
nucleosynthesis scale because baryogenesis forbids a too high entropy release during the decay 
of the scalar field (AS'<10^) in order not to wash-out any pre-existing baryon asymmetry 
0, turning, ultimately, into a further lower bound on the mass which should be larger than 
A possible solution to the Polonyi problem could seem to fine-tune the initial 
amplitude of the classical background oscillations to be much smaller than 1 (in Planck 
units). 

In this case on top of the constraints on the classical scalar background one should 
also be aware of other possible bounds coming from the scalar field fiuctuations whose non- 
relativistic modes might be amplified, for instance, during some infiationary epoch preceeding 
the coherent oscillations regime 0. 

String moduli are coupled to Standard Model fields only through Planck scale suppressed 
interactions and their effective potential is exactly flat at small coupling in the supersym- 
metric limit, but can become curved thanks to non-perturbative effects or because supersym- 
metry breaking. String moduli, then, also suffer of the Polonyi problem 0. Among the 
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moduli fields a very special role is played by the dilaton background which provides a unified 
value for the gravitational and gauge coupling constant at the String scale A5 and which 
then couples not only to the Einstein-Hilbert term but also, for example, to the abelian and 
non-abelian gauge fields ||^ present in the theory. Precision tests on the equivalence principle 
would seem to imply a lower bound for the dilaton mass m (m > 10^"^ ev) and, unless we 
do adopt the approach of |]lOl (where a very fight dilaton mass is postulated), the problem 
of the coherent oscillations persists. 

One possibility would be to fine-tune the initial amplitude of the dilaton to be much 
smaller than one (in Planck units) and then the possible constraints coming from the am- 
plified dilaton fiuctuations seem quite alleviated provided they are amplified with spectrum 
increasing in frequency as it seems reasonable in a string cosmological scenario [0. 
There have been recently few interesting suggestions in order to solve the moduli problem 
. A second stage of infiation able to dilute the energy of the moduli field was suggested 
13| and, more recently, it was also pointed out that a mass term contribution of the order 



Cm?H'^{t) (with C >> 1) could smear the amplitude of the initial coherent oscillations 
definitely relaxing the problem [^]. 

Our purpose in this paper is not so ambitious and we will not specifically address the 
moduli problem. We want only to point out that the amplification of the gauge fields 
excited by the coherent oscillations of the dilaton coupling should also be taken into account 
in the same way as the scalar field fiuctuations were shown to be important in the 
standard Polonyi problem. It is actually well known that in General Relativity the only 
fields which can be amplified through the evolution of time- dependent background geometry 
are the ones whose equations of motion are not invariant under the Weyl rescaling of the 
metric tensor (such as, for example, scalar and tensor inhomogeneities [l^) and since the 
evolution equations of the Maxwell fields (and of the chiral fermions) are invariant under the 
conformal rescaling of the metric background, it was also correctly concluded that neither 
the electromagnetic inhomogeneities nor the chiral fermions can be successfully amplified. In 
string inspired backgrounds the situation changes and gauge fields fiuctuations might grow 
thanks to the dilaton coupling in the same way as the coupling to the background curvature 
allows the parametric amplification of gravity waves in General Relativity. 
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The amplification of tlie electromagnetic fluctuations was discussed in the context of 



the magnetic field generation []T3|, but not studied, so far, in the context of the coherent 
oscillations of the dilaton field. 

In the low energy string theory effective action the evolution equations describing 
the dynamics of the massless modes can be written either in the Einstein frame (where the 
Einstein-Hilbert term and the dilaton field are formally decoupled and the Planck scale is 
really a constant) or in the String frame (where the coupling of the dilaton to the Einstein 
Hilbert term is explicit and the typical scale of the theory is the String scale [while the 
Planck scale is no longer a constants since it evolves in time as Ap = ea A5]). As long as 
the corrections in the string tension are not included the two frames are equivalent (up to a 
conformal transformation and to a ri-definition of the dilaton field). Moreover it was checked 
that the amplitude of the excited inhomogeneities is exactly the same in both the description 
even though the treatment in the String frame can be, in principle, more cumbersome [p^ . 
Since we are mainly concerned, in the present discussion, with small curvature scales [H < 
Mp) and small dilaton couplings (e*^/^ <^ l)we find useful to work, from the very beginning, 
in the Einstein frame, where the low energy effective action takes the form: 

S = -J d'x^{R - \g^'d,<pd,<P + V{<P) + \e-^F;,F'^>^'' + ...) (1.1) 

(we adopted "natural" units by setting IGvrG = 1 and we neglected the possible contribution 
of the antisymmetric tensor field which is not directly relevant for the problem which we are 
discussing; is a generic gauge field which will be considered abelian; = $ — In Vg = In (7^ 
is connected with the four dimensional gauge coupling [$ is the 10- dimensional dilaton field 
and Vg is the volume of the 6-dimensional compact internal space] ). 

The origin of the dilaton potential V^(0) is non perturbative and related to supersymmetry 
breaking. On theoretical ground, in critical superstring theory the dilaton potential has to go 
zero as a double exponential as — —00 (weak coupling), namely V(0) ~ exp (— exp (—0)) 
with a positive (model dependent) constant. On more physical ground V{(j)) is beheved 
to have also one (or more) minima for some finite value of the dilaton ~ 0mm- In general 
we can say that after a phase of decreasing coupling driven either by the dilaton's kinetic 
energy or by a radiation dominated background the coherent oscillations start at a curvature 
scale if ~ m ITT], |12|. 
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The fluctuations of tlie gauge fields can then be amphfied in two different ways : either 
through the oscillating phase itself (for Fourier modes k<m) or through the sudden change 
from the (model-dependent) pre-oscillating epoch to the oscillating one. In the first case 
the gauge fields are amplified in a resonant way with a mechanism very similar to the ones 
recently introduced in the studies of the so-called explosive infiaton decay [l^, |l^] during the 
pre-heating dynamics |]16|, |l^, |l^, In the second case the gauge fields are amplified in a 
non-resonant way since the main effect will be given not by the oscillations themselves but 
by the the change of the background evolution towards a phase effectively dominated by a 
non-relativistic pressure-less fluid which ends up, ultimately, with the final dilaton decay at 
a curvature scale ~ F ~ m^/Mp. The amplified electromagnetic inhomogeneities should 
not over-close the universe and then the energy density of the fluctuations should be always 
smaller than the energy density of the relativistic plasma. The critical energy density bound 
applied on the amplified fluctuations can then be translated into a bound on the parameters 
of the background during the phase of the coherent oscillations, leading to possible bounds 
on the dilaton mass. 

The outlined picture of the dilaton relaxation might be unreliable for many reasons. By 
studying the coherent oscillations only in the vicinity of the minimum of the potential we 
will tacitely assume that the potential is only quadratic. Detailed studies of the preheating 
dynamics both in flat |^ and curved space showed that this approximation is quite 
misleading. The quadratic potential is in fact, by itself, unstable under renormalization and 
a is needed from the one-loop divergences. In the case of the present paper, moreover, 
even higher powers are expected given the double exponential form of V{(f)). Our results 
will not have the same degree of accuracy of Ref. ^0, 21] also because we do not include 
the back-reaction effect of the produced fluctuations on the dynamics of the background as 
carefully done in ||21| in the case of de Sitter space. 

The case which we are going to discuss is also quite peculiar because of the exponential 
coupling of the dilaton to the gauge fields which was not specifically considered in any of 
the previous discussions of the preheating dynamics. This limitations will be reflected in our 
calculations and we will come back to this point later. 

The plan of this paper is then the following. Sec. II is devoted to the study of the resonant 
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amplification of electromagnetic fluctuations occurring for modes k>m, while in Sec. Ill we 
will turn our attention to the non-resonant case {k<m) and we will show, within some very 
elementary example, how the bounds on the amplified electromagnetic fluctuations can be 
turned into bounds on the dilaton decay scale and then on the dilaton mass. 



2 Resonant Amplification 

In a conformally flat metric of FRW type 

g^, = dtag{l, -a\t), -a\t), -a\t)) (2.1) 

the background equations derived from the action ( |1 . 1| ) become 

= ^{^' + 2V), H = -\^' 
dV 

<i) + 3/70 + — = (2.2) 

{H = (In a)', ■ = d/dt). For scales H<Hmax < Mp ( where Hmax is the maximal curvature 
scale during the pre-oscillating phase ) the evolution of the dilaton can be quite different 
depending upon the particular model (see for instance 0, |Tl|, [12|). Our discussion does not 



rely on any special background during the pre-oscillating phase but only on the existence of 
such a phase, occurring for H > m. We will only assume that the curvature scale and the 
coupling constant {g = e^) are both decreasing in time (if < 0, ^ < 0). 

For example in the post-big-bang evolution of the String cosmological models the back- 
ground (right after Hmax) is dominated by radiation {H ^ l/2t) and the dilaton decreases 
according to the third equation in ( |2.2| ) reaching a constant value before the onset of the 
coherently oscillating regime, i. e. = + (j)i{H / HmaxY^"^ ■ If the background is dominated, 
during the pre-oscillating phase, by the dilaton's kinetic energy 0, |12| {H ~ l/3t, p<^ ~ a~^) 
the dilaton decreases as = 0o + 0i In {Hmax/ H). The potential appearing in Eq. (|2.2|) can 
be approximated (according to the considerations of the Introduction) by 

2 

V{<P) = —{<P-<Pminfe^V{-C^^^V{-<t>)) ■ (2.3) 

In the following discussion we will keep only the quadratic terms in the potential. This 
approximation will allow our estimates but at the same time we stress again that already 
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the (flat space) solutions of the theory including quartic self-interactions might turn out 
qualitatetively and quantitatively different giving rise to sizeable errors (up to 40%) in the 
number of produced particles |^ . 



As soon as H ^ 0{m) the dilaton starts oscillating with initial amplitude 0o — 0{1) in 
Planck units. Moreover, for interactions of gravitational strength we also expect (pmin — 0{1) 
and, without fine-tuning, |0o — 0mm| — Oil). 

In the oscillating regime a self-consistent solution of the Eq. ( ^.21 ) with the potential 
(PD is: 

,,, , 4 .(hn cos mt + (hn smmt (ht sin mt cos^ mt + 6 ^ cos mtsiv? mt. ^,,H.n. 

m = +^{- + r-i^ + c - 

V3 rnt [mty m 

(2.4) 

"W-"»("''^'^'l^^i^-2^"°«^''>' (2.5) 
(with |0+|2 + |0-|2 = a2). 

The leading term (in H/ m) of Eq. ( |2.5| ) turns out to be more important in the case of the 
non-resonant amplification which rely on the change in the background evolution among the 
pre-oscillating and the oscillating phase. It is easy to transform the previous expressions to 
conformal time 7] {a{ri)dri = dt) by keeping only the leading order behaviour in H/m in Eq. 
( p. 51 ). We will use the conformal time in order to describe the non-resonant amplification 
(where the expansion of the universe is mainly concerned). On the other hand it is far 
more convenient to analyse the resonant amplification in cosmic time where the solutions of 
the background equations can be directly inserted into the gauge fields evolution equations, 
leading, finally, to more tractable expressions. 

From now on we will specialise our discussion to the case of abelian fields ( but our 
discussion can be easily extended to the non-abelian case provided e*^/^ <^ 1). The evolution 



equation for the Fourier modes of the abelian gauge fields (derived from the action (|1.1|) and 
in the gauge Aq = 0, VjA* = 0) are, for each of the two physical polarizations 

Mk)" + [k' - g{g-r]Mk) = 0, g{4>) = J (2.6) 

(where Ai = g{4>)Ai and ' = d/drj). It is useful to rewrite the same equation in cosmic time 

Mu;) + [iu'-Git)]f,iiu) = (2.7) 
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where uj = k/a, fi = \faAi and 



(the second exact equahty follows directly by repeated use of the background equations 
( p.2| )). By inserting Eq. ( |2.4| ) in Eq. ( p.8| ) and by keeping only the leading terms in H/m 
we obtain that Eq. (|2.7| ) can be rewritten as 

^ + ['^e//(^)-2ee//cos2^]/, = 0, z = — (2.9) 

where Sefj = 4(^)^(2;/2)~''/^ , e^fj = ez~^ (e = 8/(3-\/3)0d ; in order to write the previous 
equation we have chosen, for simplicity, that 0o = in Eq. ( p.4|) and we also assumed 
(pmin — 4'o)- If we focus our attention on the regime where k » m (i. e. H » m) we could 
neglect, in the first approximation the universe expansion and then 6eff ~ 5 = 4(A;/m)22^/^ 
eg// — e. In this resonant approximation Eq. ( |2.9| ) becomes exactly the well known Mathieu 
equation ||26|. Notice that when e = and 6 = 2n + l (ra = 0, 1, 2...) there is an exact solution 
of Eq. ( |2.9|) with period 27r. U 6 = 2n there is an exact solution with period vr. We can then 
say that the instability boundaries in the 6 — e plane will cross e = at the points 6 = {2n)'^ 
and S{2n+ 1)^ which are the conditions for the parametric resonance to hold. This analysis 
can be pursued very easily if e <^ 1 (i.e. 0o ^ 1). 

A more precise map of the boundaries of the unstable region can be obtained in the 
6 — e plane to the desired order in e. The analytic treatment shows that in this case the 
solutions of Eq. ( |2.9| ) are exponentially unstable but only in narrow resonance bands (see, 
for instance, [0 |T9l and |25|) \E„] - n2"23-2"[(n- l)!]-ie"5i-" with 6 = riT ^ and En < 6. 
Each resonance band with n > 1 corresponds to keeping terms 0{e'"') and then the width of 
each successive resonance shrinks but near the center of each resonance band the solutions 
take the form fi"'\z) = (45(A;))~-'^/^ exp {±a^^^z) (with a^^Veal and positive). 

As it was extensively noticed in the context of the studies of the preheating dynamics in 
the narrow resonance regime the amplification of the inhomogeneities can only be very mild, 
bmce m our case e = 8(3V3)-i ~ 1.53 we might argue that, provided we do not fine-tune 
(pQ <^ 1 we are not in the narrow resonance regime but in the broad resonance. We notice 
that the mode function G{t) appearing in Eq. ( |2.9| ) reduced simply to a damped cosinus only 
because we kept only quadratic terms in the dilaton potential in the vicinity of its minimum. 
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As shown in already the inclusion of term (required by stability under renor- 
malisation) leads to a Lame' equation where the cosinus in the mode function is replaced 
by an elliptic cosinus. Needless to say that the most accurate analysis is performed using 
the Lame' equation. In the present case, however, the mode function would look even more 
complicated than the one of the usual Lame' equation (even if we keep only quartic self- 
interactions in V{(j))), because of the exponential coupling of the dilaton to the gauge fields 
which is very different from the usual couplings discussed in the context of the preheating 
dynamics [|16|, [l^, |g. 



The previous (and following) estimates based on Eq. ( |2.9| ) are then expected to be 
drastically changed in a fully consistent non-linear calculation and we only use it as a naive 
toy model but, none the less we find useful to illustrate with simple examples the effect we 
are discussing. 

If in the case e ~ 1 we can have a resonance (in flat space) this can be expected to 
be shut off pretty soon by the universe expansion. A simple Runge-Kutta method can be 
employed after the second order differential equation (|2.9| ) has been reduced to a system of 
two first order in the variables Q = fi and P = We took (pQ = 1 and e = 8(3-\/3)~^ and 
we explored different values of k/m after having normalised Q at the value of the quantum 
mechanical fluctuations {Q{zo) = /{zq) = {ASeff{zQ))~^^^), and Zq = 2 [i.e. to = '^/m] is the 
initial integration time). The results illustrated in Fig. 1 emphasize the importance of the 
universe expansion. We plot fi{z) for different values of the parameters k/m and A = (pQ . 

The growth of the solutions in fiat space and the consequent particle production is seri- 
ously affected by the universe expansion. The solutions discussed in this section neglected 
completely the problem of the back-reaction of the produced particles on the evolution of the 



geometry. It was actually stressed in ||I6|, that the effect of the produced inhomogeneities 
modifies the effective equation of state from p = (in the regime of coherent oscillations) 
to p = p/3 (when the produced fluctuations dominate). In our case we can estimate that 
the back-reaction effect can be safely neglected for times At ~ l/m(A + 0{H/m)) where A 
is a number which depends upon the particular model of dilaton relaxation and upon the 
amplitude of the dilaton (in Planck units) at the beginning of the oscilating phase. 
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Figure 1: We plot few solutions of Eq. (2.11) for A = 0q = 1 (i.e. e = 1.53) in the resonant 
case {k > m). The effect of the universe expansion on the fiat space solutions is described. 

3 Non-resonant amplification 

We just showed that the amplification due to parametric resonance seems quite small and we 
now turn to the problem of non-resonant amplification occurring for modes k<m. We will 
discuss in this section a "minimal" toy model where the dilaton rolls down during a (model 
dependent) pre-oscillating phase, settles to its minimum and then decays at a curvature scale 
determined by its decay width Hd ~ F ~ m? /Mp producing a radiation dominated epoch. 
We are then led to study the amplification of massless gauge fields in the following model : 

g = = (^)/^, a{r^) ~ (^)^ v < Vd, and a(r/) ~ (^), > r/, (3.1) 

^max ^max 

(with (3 <0 and a < 0). 

According to this simple model of background evolution the effective potential V{r]) = 
g{g^^)" appearing in Eq. ( p.6|) will be : 

^ = ^(^±11, ^ < and V{r]) = 0, r/ > r/^ (3.2) 

A particular solution of Eq. (|2.6| ) for Ai{k) can thus be written in terms of the first and 
second kind Hankel functions if'-^-* and H^"^^ corresponding to free oscillating modes in the 
|?7| oo limit. 
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The effective potential barrier leads to an amplification of the initial fluctuations, ore 
equivalently, to particle production from the vacuum. Starting with an incoming Fourier 
mode which is of positive frequency with respect to the asymptotic state at the left of the 
barrier (77 —00), one has in general, for — 00 a linear combination of modes which are of 
positive and negative frequency with respect to the new vacuum. The mixing coefficients c± 
define the Bogoliubov transformation connecting the two vacua and determine the spectral 
distribution of the amplified fluctuations. 

The general solution of Eq. ( p.6|) in the two distinct epochs defined in Eq. (|3.1| ) - (|3l^ ) 
can be written for each Fourier mode Ai{k) 

Mk) = ^^Hj^\kr]), ri<r],, and = ^(c+e"^'^'^ + c.e*'^") r] > r], (3.3) 

The Bogoliubov coefficients c± can be fixed by the two conditions given by matching Ai and 
A'j^ in 1] = rjd. It is well known that the Bogoliubov coefficients determined in this "sudden" 
approximation can, in principle, lead to ultraviolet divergences in the energy density of 
the amplified fluctuations. The physical reason for this is that the modes with frequency 
k > l/rjd the sudden approximation cannot be rigourously applied and the calculation should 
be performed with a smooth effective potential barrier interpolating among the two vacua. 

To smear the potential step will produce an exponential suppression of |c_| for modes 
k > l/rjfi. Then in the frequency band we are interested in the amplification coefficient will 
be then given by |c_| ~ (/cr^^)^^^''^' and the energy density of the produced massless gauge 
bosons p, per logarithmic interval of proper frequency cilna;, is obtained by summing over 
the two polarisation states of the amplified waves: 

dp 



dhiuj 



2Nujy_\^ (3.4) 



We included also the factor which counts the different massless gauge fields appearing 
originally in the low energy string theory effective action (|1 . 1| ) which can be, in principle, 
also quite large, of the order of 10^ — 10^ (as it happens, for example, in heterotic string 
theory [^). In order to obtain the final expression of the spectral energy density of the 
amplified fluctuations in critical units we divide Eq. (|3^) by the radiation energy density 
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produced by the dilaton decay (pc = H'j^^pi.'^d/ o^Y) and we obtain 

r]M = l-^^iv(^)-%i(^)^ (3.5) 

Pcdmuj ujd Mp 

The spectral energy density of the amphfied inhomogeneities has to be always VL{uj,t) < 1 
for all the frequencies and for all the times. If the energy spectrum is either flat or decreasing 
{\^A — 3/2) the critical density bound should be applied at the lowest amplified frequency. 
If the spectrum is increasing (|yu| < 3/2) the bound should be applied to the highest possible 
frequency which coincides, in our case with uj ~ uj^. Therefore the critical density bound 
implies, respectively, 

^.<iV-|(^)^, for > ^ and, ^<N~K for |/.| < ^ (3.6) 



Mp 'ud - 2 ' 

where uji is the infrared scale of the spectrum. As an example we consider the case where 
during the pre-oscillating phase the evolution is driven by the kinetic energy of the dilaton 
starting from the scale Hmax < Mp. In this case /? = —VS and \n\ = (-\/3 — l)/2. The 
critical density bound applied to the highest amplified frequency gives from ( |3.6| ) in terms 
of the dilaton mass m<N~^^^Mp (we used Hd — {m? /Mp). It is amusing that the critical 
density bound implies a constraint for the dilaton mass which depends upon the number of 
the amplified massless gauge fields. If ^ 2 x 10'^ this bound basically implies m<10^^Mp. 
Our example shows clearly that, in spite of the simplified model of background evolution it is 
possible to use the non-resonant amplification of massless gauge fields in order to constrain 
the dilaton mass. Since our purpose, in this section, is only illustrative, we only mention 
that it is possible to obtain more complicated exclusion plots either by leaving /x as a free 
parameter, or by adding to our "minimal" example the effect of the intermediate phase 
(occurring for Hd < H < m). We are anyway studying modes k < m and the addition of 
this phase will only change the Bogoliubov coefficient |c__| (and then the energy density of 
the amplified fiuctuations) through the inclusion of an extra parameter (the duration of the 
intermediate phase) which will make broader (in conformal time units ) the "under-barrier" 
region in the effective potential appearing in Eq. ( |3.2| ) . In these cases further constraints 
relating the dilaton decay scale to the other parameters of the background model can be 
obtained. 
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We discussed the resonant and non- resonant amplification of massless gauge fields during 
a phase of dilaton oscillations and our conclusions are that the resonant amplification is a 
less efficient mechanism unless the value of the dilaton at the beginning of the oscillating 
phase is fine-tuned to be orders of magnitude bigger than the Planck mass which appears to 
be unreasonable for interactions of gravitational strength. The important caveat is, anyway, 
that the estimates we presented were done in some oversimplified model of dilaton relaxation 
and, for a more refiable answer, a complete non-linear analysis would be needed. The non- 
resonant amplification of massless gauge fields seem, already in our example, a more efficient 
mechanism and can very likely lead to possible constraints on the value of dilaton decay 
scale and on the value of the dilaton mass. 
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